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ABSTRACT: We first sought to clarify the effects of discounted rate, survival rate, and lactation persistency as a component trait of 
the selection index on net merit, defined as the first five lactation milks and herd life (HL) weighted by 1 and 0.389 (currently used in 
Japan), respectively, in units of genetic standard deviation. Survival rate increased the relative economic importance of later lactation 
traits and the first five lactation milk yields during the first 120 months from the start of the breeding scheme. In contrast, reliabilities of 
the estimated breeding value (EBV) in later lactation traits are lower than those of earlier lactation traits. We then sought to clarify the 
effects of applying single nucleotide polymorphism (SNP) on net merit to improve the reUability of EBV of later lactation traits to 
maximize their increased economic importance due to increase in survival rate. Net merit, selection accuracy, and HL increased by 
adding lactation persistency to the selection index whose component traits were only milk yields. Lactation persistency of the second 
and (especially) third parities contributed to increasing HL while maintaining the first five lactation milk yields compared with the 
selection index whose only component traits were milk yields. A selection index comprising the first three lactation milk yields and 
persistency accounted for 99.4% of net merit derived from a selection index whose components were identical to those for net merit. We 
consider that the selection index comprising the first three lactation milk yields and persistency is a practical method for increasing 
lifetime milk yield in the absence of data regarding HL. Applying SNP to the second- and third-lactation traits and HL increased net 
merit and HL by maximizing the increased economic importance of later lactation traits, reducing the effect of first -lactation milk yield 
on HL (genetic correlation (iq) = -0.006), and by augmenting the effects of the second- and third-lactation milk yields on HL (Tq = 0.118 
and 0.257, respectively). (Key Words: Economic Weights, Milk Yield, Persistency, Herd Life, Selection Index) 



INTRODUCTION 

High milk production and long herd life (HL) are 
important to the dairy industry. In Japan, the average 
number of lactations of a commercial dairy cow is about 3. 
Lifetime milk yield, in principle, is defined as a total milk 
yield until the time the cow is culled but can be defined as a 
total milk yield from as many as five lactations. Previous 
selection index coefficients of the first three lactation milk 
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yields to maximize the first five lactation miUc yields 
(Weller et al., 1984) ignored the differences in reliability of 
estimated breeding value (EBV) on various traits of the first 
three lactations. However, selection index coefficients can 
be affected by the reliability with which EBV predicts a 
particular component trait (Togashi et al., 2011). The 
reliability of EBV differs with selection traits of the first 
three lactation milk yields and lactation persistency as well 
as between bulls and dams selected as breeding candidates 
(Hagiya et al., 2010). Genetic superiority of animals 
selected as breeding candidates or genetic difference 
between selected and all the animals is estimated more 
precisely by allowing for differences in the reliabilities of 
index component traits. 

Prolonging the HL of high-producing cows decreases 
replacement and health costs and thus increases profit. Herd 
hfe is a complex trait that is determined by survival rate, 
which in turn is influenced strongly by multiple factors such 
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as milk production, disease, conformation, selection 
intensity, market conditions, and so forth. Obviously, 
improving the survival rate would increase HL. Multi-trait 
selection of breeding candidates in which HL is included as 
one of the component traits requires knowing the economic 
weight of HL in relation to that of all other component traits. 
Estimating the relative economic weight of HL requires 
milk, carcass, and calf revenues and feed, health, labor, and 
opportunity costs (Van Arendonk, 1991; Van Raden and 
Wiggans, 1995; Weigel et al., 1995; Jagannatha et al., 1998). 
Currently, the Japanese scheme of selection for lifetime 
performance is based on a linear combination of lifetime 
milk yield and durability (an indicator of HL), in which 
these traits are weighted by 1 and 0.389, respectively 
(NLBC and LIAJ, 2010). 

Therefore, net merit is defined as a linear combination 
of total milk yield from the first five lactations and HL, 
weighted by 1 and 0.389, respectively, in units of genetic 
standard deviation. Thus, comparison of selection index is 
investigated based on the magnitude of net merit. 

Improving lactation persistency improves animals' 
utiUzation of cheap roughage (Solkner and Fuchs, 1987), 
reduces stress due to high peak yield (Zimmermann and 
Sommer, 1973; Jakobsen, 2000; Yamazaki et al., 2010), 
increases profit (Dekkers et al., 1998), and prolongs HL 
(Weller et al., 2006). We therefore evaluated the effects of 
lactation persistency and milk yield during the first three 
lactations as component traits of the selection index on net 
merit. 

Discounted rate is the quantification of the response to 
selection according to the discounted cash-flow procedure 
(Hill, 1971). In the current study, we investigated the effects 
of discounted and survival rates on net merit. We 
hypothesized that survival rate would increase relative to 
the economic importance of later compared with earher 
lactation traits and would increase the accumulated 
economic value of the first five lactation milk yields during 
a given period. However, the reliability of EBV is lower for 
later compared with earlier lactation traits (Hagiya et al., 
2010), making it essential to increase the reliabiUty of EBV 
for later lactation traits to make good use of increased 
economic importance in later lactation traits due to increase 
in survival rate. Single-nucleotide polymorphism (SNP) has 
been applied to increase the reliabiUty of EBV in young 
bulls (Van Raden et al., 2009), and the genetic correlation 
between HL and first-, second-, and third-lactation milk 
yield is -0.006, 0.118, and 0.257, respectively (Hagiya et al., 
2010). We therefore assessed whether application of SNP to 
improve the reliabilities of EBVs for HL, lactation 
persistency and milk yield in later lactations consequently 
increased net merit and HL, reduced the effect of first- 
lactation milk yield on HL, and augmented the effects of 



second- and third-lactation milk yields on HL. 

MATERIALS AND METHODS 

Derivation of survival matrix 

We divided the life cycle of a dairy cow into 5 stages: 
stage 1 refers to the period from birth to the end of the first 
lactation; stage 2 is from the end of the first lactation to the 
end of the second lactation; and so forth. The survival 
matrix Sij then is derived, in which the row (i) represents 
the lactation number and the column (j) represents the age 
(in months) at the end of the i* lactation. Thus, matrix S 
represents the age distribution of a cow population. 
Currently in Japan, the average age at first calving is 25 
months, and the calving interval is 14 months (MAFF, 
2011). Let roi be the survival rate from birth to the end of 
first lactation or the end of stage 1. Let be the survival 
rate from the end of the i* lactation or the end of the i* 
stage to the end of the j* lactation or the end of the j* stage, 
e.g., ri2= survival rate from the end of the first lactation or 
the end of the first stage to the end of the second lactation 
or the end of the second stage. The survival rate (ry) is 
determined by a cow being alive and getting pregnant from 
the end of the i"" lactation to the end of the j* lactation. The 
survival rates of roi, rn, '"23, ''34, and currently in Japan are 
0.6900, 0.7773, 0.7525, 0.5936, and 0.5464, respectively 
(Hagiya et al., 2010). The survival matrix starts when a 
female calf is born. The investment period is set at 120 
months. The elements of Si j (i = i_ 2, 3, 4, 5 staged = 1,2, 120 months) 
are derived as described by Weller et al. (1984). Increased 
survival rates are set as t times as large as the current 
survival rate (t = 1.05, 1.1, 1.15, 1.2, and no loss rate of 
100%). 

The relative economic weights of the first five lactation 
milk yields 

Genetic Ufetime milk yield (H5) is defined as follows: 
Hj = WjGi + W2G2 + W3G3 + W4G4 + W5G5 

where Gi and Wi are the genetic value and relative 
economic weight of the i* lactation milk yield, respectively 
ii = 1 to 5). H5 is genetic Ufetime milk yield derived by 
using only milk yield for the first five lactations. 

The discounted expressions (Hill, 1971; Weller et al., 
1994) are as follows: 



' (1 + 6)^ 

where d, is the appropriate discounting factor, for which the 
subscript j is the time interval (in months) from the original 
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investment to mean trait expression (j = 1, ... , 120 months), 
and 5 is the discounted rate per month %, e.g., 5 = 0/12%, 
1/12%, or 4/12% per month. The relative economic weights 
for the first five lactation milk yields (w'5 = (wi, W2, W3, W4, 
W5)) reflect the actual average of each lactation milk yield, 
so that W5 is Wj = AMj^jSj^j^^d with the vector d = (di, ... 

di2o)' and the diagonal matrix AM5X5 whose diagonal 
elements are the averages of the i' lactation milk ( M; , ' = 
1 to 5). The discounted total milk yield during the given 
period of 120 months from the original investment is 

described as Wj 1 . 

Net merit defined in terms of the first five lactation milk 
yields and HL 

The relative economic weights (in genetic standard 
deviations) of total milk yield (the sum of the first five 
lactation milk yields) and HL were studied as: 

wl„k ■ w'hl = 1 : 0.389 , where and W^^ are the 

relative economic weights (in genetic standard deviations) 
for the first five lactation milk yields and HL, respectively. 
The ratio of w^-,^ : = 1 : 0.389 is equivalent to the 

ratio of 7.2:2.8 currently used in Japan and is in close 
agreement with the reported value of 7.0:2.5 (Burnside et al., 
1984; Allaire and Gibson, 1992; Weller and Ron, 1992; 

Dekkers, 1994). Net merit {^„et merit )> defined in terms of 

the first five lactation milk yields and HL in genetic 
standard deviations, is as follows: 

H* 

net merit 

_ o-H3(WjGi +W2G2 + W3G3 + W4G4 +w;G5) 

where w?* {i = 1 to 5) is the relative economic weight for 
the i* lactation milk yield described as vv'" =w, /(^W;); 

i=l 

5 

^ vf ° = 1 > wl is the relative economic weight for HL, 

i=i 

cTh5 = (Variance(wJGi + w^G^ + w^G, + w^G^ + w^Gj))"^ 
and cr^ is the genetic standard deviation of HL. On the 

<-'HL ° 

Other hand, / w* ■„ = = — — — ^ . so that net merit 
7-2 cr^5 



(H ngtmerit^ iu uuits of actual values of milk (kg) and HL 
(days) is as follows: 



= w"G^ + w^G, + w^Gj + w^G^ + W5 G5 + WgGaj, 

Comparison of selection indices in terms of maximizing 
net merit 

Lactation milk yield and persistency during the first one, 
two, or three parities are taken as the component traits of 
selection indices (li-L) to maximize net merit (H ) 

and to clarify the effects of adding lactation persistency as a 
component trait of selection index to selection index 
consisting of only milk yields. Selection index (I7), 
comprising the same components as those for net merit, is 
set as a basis of comparison: 

l3=biM,+bjM2^ 

14 = biMj + bjMj + bjPi + b^Pj 

15 = bjMj + bjMj + b3M3 

I, = biMj + b,M, + b3M3 + b^Pi + b^P, + b,P3 ^ 
I7 = biMj H-bjMj +b3M3 +b4M4+b5M5 + b^hl 

where M;, Pj , and hi are the EBV for the i* lactation milk 
yield, i* lactation persistency, and HL, respectively and bi is 
the selection index coefficient for the i* component trait of 
selection index. Lactation persistency is defined as the 
difference in EBVs for 240 DIM and 60 DIM (Togashi et al., 
2008). The selection index is derived as described by 
Togashi et al. (2011), in which differences between the 
reliabilities of the EBV of the component traits of the 
selection index are taken into consideration. 

Here we use the index Ig = b|M|+b2M2+b3M3+b4Pi+b5P2 
H-bgPs as an example to demonstrate the derivation. In 
matrix notation, the index can be written as, 

l6= h\g„ 
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where b, = (bib,b3b,b,b,) and g, = (M,M,M3P,P,P3). 

As indicated earlier, the goal is to maximize the genetic 
superiority of animals selected as breeding candidates in 
terms of net merit {H . ). The covariance matrix 



between Ig and H 



is as follows: 



where w' = {w" w" w" w" w" w" ) and 



'gi'g4'^G1 
^GP'J^GA^W 



'G2'G5^GXaa. 

''Gi'GS^a^S 
'"(S'l^CS^dil.GS 



'CP rG2" (3*1,02 ' CpI'G3^CP1,G3 'GPl'GA^&\.Gi 

^Ca'2^G2^a'2,G2 ^GP2^Gi^GP2,G3 'ffl'2'c4^GP2,G+ '&2'G5^GP2.G5 



't3'3'G2"t3'3.G2 'GP3'G3^GP3,G3 



'GP3'G5"ra'3.G5 



'G2'(aiL^G2,GM. 
^GS^Caa.'^GS.GHL 



'GP2'oiL^GP2,GliL 
^GP3^0IL^CPi,CHL 



The variance of Ig is as foUows: 



al =b:P.b. 



where GP, is the genetic value of the i* lactation persistency, 
is the accuracy of the i* lactation milk yield EBV, and 

r^yp. is the accuracy of the i* lactation persistency EBV. 

The reliabilities of the EBV of selection traits of candidate 
bulls and the genetic parameters were derived from Hagiya 
et al. (2010) and Togashi et al. (2008). The square of the 
accuracy is rehabihty. The average reliability is used in this 
study. 

The genetic superiority of selected bulls as breeding 
candidates or the genetic difference between candidate bulls 
and all the bulls due to selection on Ig for net merit 

(^„«™.w<)is as follows: 



where Oj^ is the standard deviation of the index Ig and i 

is the selection intensity. 

And selection accuracy is as follows: 



where G is the genetic covariance matrix of the first five 



lactation milk yields and HL. The genetic superiorities of 
selected buUs as breeding candidates in the first five 
lactation milk yields (AGi,AG2,...,AG5) and HL (AGj^^) 
are shown as follows: 

[AG, AG, AG3 AG, AG3 AG„j'=C,bJ/c7,,. 

We compared the seven selection indices in terms of net 
merit, selection accuracy, genetic superiorities over the total 
of the first through fifth lactation milk yields (^ag. 

herd life ( A//L ). Selection intensity was set to unity {i = 1 ) 
for all selection criteria compared in the current study. 

Application of SNP 

Currently in Japan, the reliabilities of EBV among 
selected bulls as breeding candidates for milk yield of the 
first five lactations are 0.816, 0.708, 0.597, 0.550, and 
0.500, respectively; those for persistency of the first three 
lactations are 0.689, 0.612, and 0.493, respectively; and that 
for HL is 0.570 (Hagiya et al., 2010). The reliabilities of 
first-lactation traits are higher than those of second or third 
lactation traits, because collection of records is more rapid 
for first-lactation compared with second- and third-lactation 
traits. Using SNP to improve reliability is more important 
for second- or third-lactation traits than first-lactation traits, 
due to this difference in data collection time. Therefore, we 
considered three scenarios: increasing the reliabilities of 
second- and third-lactation traits and of HL by 10% while 
maintaining the reliabilities of first-lactation traits at their 
current levels; increasing the reliabilities of second- and 
third-lactation traits and of HL by 20% while maintaining 
the reUabiUties of first-lactation traits at their current levels, 
and increasing the reliabilities of second- and third-lactation 
traits and of HL by 20% and the reliabilities of first- 
lactation traits by 10%. 

RESULTS AND DISCUSSION 

Effects of survival and discounted rates on tlie relative 
economic weights of the first five lactation milk yields 

The relative economic weights for the first five lactation 
nmlk yields in terms of actual milk yield (kg) are shown in 
Table 1, which is based on survival matrix S,^- and shows 

the effects of both survival and discounted rates on the 
relative economic weights of the first five lactation milk 
yields, the pooled weights of the third, fourth, and fifth 
lactation milk yields, and these economic weights under a 
100% survival rate. The relative economic weights of the 
second, third, fourth, and fifth lactation milk yields (W2, W3, 
W4, and W5) were expressed as ratios of the first lactation 
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Table 1. Relative economic weights of tlie first five lactation milk yields 



Discounted rate i lactation milk Level of survival rate (t times as large as the current survival rate) 

%/month) (M,.kg)' 1.0 1.05 1.1 1.15 1.2 Survival rate^ 



A/1 O 


1 
1 


l.UUUU 


1 AAAA 
i .UUUU 


1 AAAA 
l.UUUU 


1 AAAA 
1 .UUUU 


1 AAAA 
l.UUUU 


1 AAAA 
l.UUUU 




z 


U.o /4/ 


A Al OA 

u.yio4 


U.ytJZZ 


1 AA^A 

1 .uujy 


1 l\A A^ 

i.U4yo 


1 1 T<Q 
I.IZJ j 




i 


U.oo4/ 


U. /j4j 


A QO'7Q 

U.oz/o 


A C\f\AQ 

U.yU48 


u.yojz 


i.ioy / 




A 

4 




A AAfYi 

U.44UJ 


A <A/^0 

U.jUOz 


U.J /o4 


U.OJ /z 


i.uyjj 




J 


A ^AOiC 

U.zUjO 


U./4/5 


A OA0 1 


A 1 
U.JjOl 


U.4zzl 


1 ATQ 1 

l.U/il 




"2 yl < 


l.Zuol 


1 /I /1 1A 


1 1 
i.OJZi 


l.ojyz 


z.U04j 


J.JJOl 




total miiK; 


1 AAAA 


1 1 

l.lzz/ 


l.zooU 


1 /I ATA 

1.4U/U 


i.j /uy 


0 CI QA 

z.jiy4 


1 /I o 


1 


1 AAAA 

l.UUUU 


1 AAAA 
1 .UUUU 


1 AAAA 

l.UUUU 


1 AAAA 
1 .UUUU 


1 AAAA 
1 .UUUU 


1 AAAA 
l.UUUU 




L 


A 1 
U.ODDl 


A AAA /I 

U.9Uy4 




A aa/;a 

u.yyou 


1 AQ AQ 

i.uiyj 


1.1142 




J 


A 

u.oDyo 


U. /JSZ 


A O 1 AO 

U.olUz 


A OOCC 


u.yo4z 


1.144/ 




4 


A O/CAI 


U.4Z/J 


A /I Al O 


A c/; 1 A 
U.0014 


U.OJ /o 


1 A/: '3 1 
1 .Uo J 1 




5 


0.1953 


0.2374 


0.2859 


0.3416 


0.4050 


1.0294 




3,4,5 


1.2340 


1.4029 


1.5874 


1.7884 


2.0070 


3.2372 




Total milk' 


1.0000 


1.1219 


1.2561 


1.4039 


1.5663 


2.5018 


4/12 


1 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 




2 


0.8403 


0.8823 


0.9243 


0.9664 


1.0084 


1.0811 




3 


0.6274 


0.6917 


0.7591 


0.8297 


0.9034 


1.0726 




4 


0.3371 


0.3903 


0.4487 


0.5128 


0.5826 


0.9710 




5 


0.1723 


0.2094 


0.2522 


0.3013 


0.3572 


0.9079 




3,4,5 


1.1367 


1.2913 


1.4600 


1.6437 


1.8432 


2.9515 




Total milk' 


1.0000 


1.1193 


1.2505 


1.3945 


1.5525 


2.4499 



' Ml =7756.9 Mj =9013.9 =9369.4 M^ =9387.4 =9196.4. ' Survival rate = 100%. 

' Total milk = Total discoiinted milk yield for the first five lactations during 120 months in compared with that of the current survival rate. 



milk ( Wj ) by setting w, =1 . The relative economic 
weights of through Wj increased as the survival rate 

increased from t = 1 to t = 1.2 for all discounted rates. For a 
given survival rate, the economic weights of W2 through 
W5 decreased as the discounted rate increased because the 
time to reach the second through fifth lactations is longer 
than to reach the first lactation. 

Among the 9 combinations determined at all levels of 
discounted rates (0/12%, 1/12%, and 4/12% per month) and 
survival rates (t = 1.0, 1.05, andl.l; t times as large as the 
current survival rate), the relative economic weight of the 
first lactation was greater than that of the second lactation, 
which in turn was greater than that of third lactation and so 
forth ( Wj > > W3 > W4 ). This effect occurred 

because survival rate decreases as parity increases, as 
mentioned earlier, i.e., the survival rates of roi, r^, r23, 
and are 0.6900, 0.7773, 0.7525, 0.5936, and 0.5464, 
respectively (Hagiya et al., 2010). The decrease in 
economic weight from the third to fifth parity was greater 
than that from the first to third parity because survival rate 
decreased more from the third to fifth parity than from the 
first to third parity. The economic weight of the fourth and 
fifth lactation milk yields increased with the survival rate. 
The relative economic weight of the pooled third, fourth. 



and fifth lactation milk yields was higher than that of first 
or second lactation milk yield, especially at low discounted 
rates and high survival rates. This result indicates that, 
compared with first- or second-lactation milk yield, the 
third-, fourth-, and fifth-lactation milk yields play important 
roles in the genetic improvement of lifetime performance, 
especially under low discounted and high survival rates. As 
an example, when the survival rate improved by 20% (t = 
1.2) under a discounted rate of 0, the relative economic 
weight of the pooled third-, fourth-, and fifth-lactation milk 
yields was about two times greater than that of the first 
parity. 

Currently in Japan, the first three lactation milk yields 
are given equal economic weights. This assignment 
overestimates the relative economic weight of the second 
and third lactations (wj and W3), particularly under the 

current survival rate (t = 1). Overestimation of W2 and W3 
decreased as the survival rate increased. In contrast, when 
the survival rate was achieved without loss (100%), the 
relative economic weight of the third lactation was greater 
than that of the second lactation, which in turn was greater 
than that of first lactation ( Wj >w^ > w, ) under discounted 

rates of 0/12% and 1/12% per month. This result occurs 
because the average yield of the third lactation milk is 
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greater than that of the second lactation milk, which in turn 
is greater than that of first lactation milk ( >m^ > M[ )■ 

The value of Wjl (the total discounted yields of the 

first through fifth lactation milks) during the 120 months 
from the original investment is shown in Table 1 as a ratio 
of the total discounted yields of the first through fifth 
lactation milks under the current survival rate. The total 
discounted yields of the first through fifth lactation milks 
during 120 months increased with increasing survival rate. 
For example, under a discounted rate of 0, when survival 
rate improved by 20% (t = 1.2), the total discounted yields 
of the first through fifth lactation milks increased by 1.57 
times compared with that under the current survival rate (t = 
1). Therefore, it is important to improve survival rate 
through both genetics and management. 

Comparison of selection indices in terms of maximizing 
net merit 

The selection accuracy and genetic superiorities of net 
merit, the total yield of the first through fifth lactation milks 
and HL derived by using selection indices (Ii to I7) for 
selected bulls as breeding candidates under the current 
survival rate and a discounted rate of 0 are shown in Table 2. 
Net merit and selection accuracy increased by adding 
lactation persistency as a component trait to selection index 
comprising only milk yields. That is, net merit and selection 
accuracy during the first; first and second; and first, second, 
and third parity components was larger in I2 (biM|+b2P|) 
than in Ij (Mi), in I4 (b|Mi+b2M2+b3P|+b4P2) than in I3 
(biMi+b2M2), and in l6(biMi+b2M2+b3M3+b4Pi+b5P2+b6P3) 
than in I5 (biMi+b2M2+b2M3), respectively. In addition, the 
genetic superiority of HL in first and second parity and in 
first, second, and third parity components was larger in I4 
(biMi+b2M2+b3Pi+b4P2) than in I3 (biMi+b2M2) and in Ig 
(b|Mi+b2M2+b3M3+b4Pi+b5P2+b6P3) than in I5 (biM|+ 
b2M2+b2M3), respectively. In particular, the genetic 
superiority of HL was 45.4 days from Ig compared with 
39.0 days from I5. The genetic superiority of the total yield 
of the first through fifth lactation milks in the first and 
second and first, second, and third parity components were 
equivalent between I3 (biMi+b2M2) and I4 (biMi+b2M2+ 



b3Pi+b4P2) and between I5 (biMi+b2M2+b2M3) and Ig 
(biMi+b2M2+b3M3+b4Pi+b5P2+b6P3), respectively. Therefore 
lactation persistency during the second and (especially) 
third parity contributed to increasing net merit, in which HL 
was increased while the first five lactation milk yields were 
maintained at the same levels derived from the selection 
index consisting of only milk yields. 

The genetic superiority of the first five lactation milk 
yields was highest for I2 (biMi+b2Pi) among the six indices 
(Ii to If,), and followed by Ij (Mi). These results reflect the 
high genetic correlations between the first-lactation milk 
yield and second- through fifth-lactation milk yields {vq = 
0.759 to 0.974), the moderate genetic correlation between 
lactation persistency during the first parity with the first five 
lactation milk yields (vq = 0.42 to 0.53), and that the 
reliabilities of the EBVs for milk yield and lactation 
persistency rank (in descending order) as Mi_ M2, M3 and Pi, 
P2, P3, respectively. 

The order of genetic superiority of HL was Ig 
(b iM ,+b2M2+b3M3+b4P ,+b5P2+b6P(,)>l5 (b ,M i+b2M2+b2M3) 
>l4 (biM,+b2M2+b3Pi+b4P4)>l3 (biM|+b2M2)>li (M,) = I2 
(biMi+b2Pi). Iiand I2 yielded the smallest AHhl, due to the 
low genetic correlations between HL and first-lactation milk 
yield (-0.006) and persistency (0.09). Tsuruta et al. (2004) 
found that the genetic correlation between productive life 
and milk yield declined from positive to 0 over the years 
from 1979 to 1993, mainly in response to changes in 
producers' culling practices. 

Dairy producers' decisions regarding culling cows 
reflect various factors such as milk yield, reproduction 
problems, mastitis, lameness, and labor costs. All of these 
factors can change over time. Functional traits (such as 
resistance to mastitis and lameness, fertility, calving ease, 
and longevity) are defined as the characteristics of an 
animal that increase production efficiency through reduced 
input costs rather than increased output of products (Groen 
et al., 1998). High milk production during the first lactation 
decreases the efficacy of the immune system in the udder 
and leads to damage that persists throughout subsequent 
lactations (Rupp et al., 2000). The genetic superiority of HL 
due to index I5 (biMi+b2M2+b3M3) or Ig (biMi+b2M2+ 



Tiable 2. Genetic superiorities of selected bulls as breeding candidates based on selection indices compared (Ii to I7)' under the current 
survival and zero discounted rates 





Ii 


I2 


I3 


I4 


I5 


l6 


I7 


Net merit^ 


594.1 


597.8 


617.1 


619.1 


657.0 


663.5 


667.5 


Selection accuracy 


0.722 


0.726 


0.750 


0.753 


0.799 


0.806 


0.811 


1 to 5 parity milk (kg) 


2,743.9 


2,753.4 


2,698.6 


2,695.1 


2,717.9 


2,717.9 


2,722.6 


Herd life(days) 


-0.6 


0.6 


14.8 


16.4 


39.0 


45.4 


46.6 



' Ii = Ml, I2 = biMi+biPi, I3 = biMi+biMa, I4 = biMi+b2M2+b3Pi+b4P2, I5 = biMi+b2M2+b3M3, 16 = biMi+b2M2+b3M3+b4Pi+b5P2+b6P3, and I7 = 
biMi+b2M2+b3M3+b4M4+b5M5+b6hi, where Mi and Pi is i* lactation milk yield and persistency EBV, respectively, and hi and bi is EBV for herd life and 
selection index coefficient for i* component trait, respectively. 

^ Net merit is a linear combination of the first five-lactation milk yields and herd life weighted by 1 and 0.389 in unit of genetic standard deviation. 
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b3M3+b4Pi+b5P2+b6P3) was about 2.5 times greater than that 
due to I4 (biMi+b2M2+b3Pi+b4P2). These results indicate 
that milk yield during and persistency of the third lactation 
play important roles in determining the duration of HL, 
because third-lactation milk yield has a higher genetic 
correlation with HL {vq = 0.257) than does either the first- 
(ro = -0.006) or second- (xq = 0.118) lactation milk yields, 
and third-lactation persistency has a higher genetic 
correlation with HL {vq = 0.210) than does either the first- 
(xq = 0.090) or second- (xq = 0.180) lactation persistency. 
The milk yield or persistency of the lactation closest to the 
culUng stage of a cow is much more important than that far 
back from the culling stage in extending the HL. 

The component traits of selection index (I7) are identical 
to those that define net merit components. This similarity 
means that net merit and selection accuracy due to selection 
index (I7) are maximal in improving net merit {H^^^ ) . 

Net merit due to Ig (b|MiH-b2M2H-b3M3H-b4P|H-b5P2H-b5P5) was 
99.4% of that of selection index (I7). Net merit, selection 
accuracy, and HL based on Ig were the largest among 
indices (Ii to l6).The first five lactation milk yields based on 
Ig was the largest among indices (I3 to Ig). The first five 
lactation milk yields based on 1 1 or I2 were larger than those 
based on I3 through Ig However, HL based on Ii or I2 was 
nearly 0. Furthermore, collecting data for HL is time- 
consuming and costly (Smith and Quaas, 1984). Therefore, 
selection index (Ig), which includes total milk yield and 
persistency of the first three lactations, as modified by the 
trait-specific EBVs, is a practical and favorable means of 
improving lifetime milk yield in the absence of data on HL. 



Table 3. Net merit, selection accuracy, and genetic superiorities of 
selected bulls as breeding candidates for first through fifth 
lactation milk yields, and for herd life based on selection index 
(Ig)' at increased survival and discounted rates 

Discounted ^ . , 



(%/month) 



Net merit^ 


0/12 


1 


663.5 




0/12 


1.2 


651.8 




4/12 


1 


665.9 




4/12 


1.2 


654.6 


Selection accuracy 


0/12 


1 


0.806 




0/12 


1.2 


0.790 




4/12 


1 


0.810 




4/12 


1.2 


0.794 


A total of the 1st to 5th 


0/12 


1 


2,717.9 


lactation imlk yields (kg) 


0/12 


1.2 


2,699.3 




4/12 


1 


2,721.6 




4/12 


1.2 


2,704.0 


Herd Ufe 


0/12 


1 


45.4 


(AHL, d) 


0/12 


1.2 


47.9 




4/12 


1 


44.9 




4/12 


1.2 


47.3 



' l6 = biMi+b2M2+b3M3+b4P|+b5P2+b6P3, where Mi and Pi are i* lactation 
milk yield and persistency EBV. respectively, and hi and bi is EBV for 
herd life and selection index coefficient for i"' component trait, 
respectively. 



^ Net merit is a Unear combination of the first five-lactation milk yields 
and herd life weighted by 1 and 0.389 in unit of genetic standard 
deviation. 

' t times as high as the current survival rate. 



Effects of survival and discounted rates on selection 
accuracy, HL, and first through fifth lactation milk 
yields 

The effects of survival and discounted rates on selection 
accuracy, net merit, HL, and total milk yield during the first 
through fifth lactations according to selection index Ig are 
shown in Table 3. As the discounted rate increased, net 
merit and selection accuracy increased but HL decreased 
because increased discounted rate puts greater relative 
economic weight on earlier rather than later lactations 
(Table 1) and because the reliabilities of the associated 
EBVs are higher for earlier rather than later lactation milk 
yields. Net merit and selection accuracy decreased sUghtly 
and HL increased when survival rate increased because 
increased survival rate puts greater economic weight on 
later compared with earUer lactations (Table 1). 

Given the current survival rate (t = 1) and a discounted 
rate of 0% per month, the genetic superiority of selected 
bulls as breeding candidates in terms of total milk yield for 
the first through fifth lactations was 2,717.9 kg, according 
to Ig. Given the same survival rate (t = 1), when the 
discounted rate increased from 0% to 4/12% per month, the 



genetic superiority in terms of total milk yield for the first 
through fifth lactations increased to 2,721.6 kg, due to the 
increased reliabiUties of the EBVs of the earlier lactation 
milk yields. Similarly, given the same discounted rate (0% 
per month), when the survival rate increased from t = 1.0 to 
t = 1.2, the genetic superiority in terms of total milk yield 
for the first through fifth lactations decreased ftom 2,717.9 
to 2,699.3 kg, due to the decreased reliabilities of the EBVs 
and increased relative economic weights of later lactations. 

Under the current survival rate (t = 1) and a discounted 
rate of 0%, genetic superiority in terms of HL was 45.4 days, 
according to Ig. Under the same survival rate (t = 1), when 
the discounted rate increased from 0% to 4/12% per month, 
genetic superiority in terms of HL due to Ig decreased from 
45.4 to 44.9 days (only approximately 1.0%). When the 
survival rate increased by 20% (t = 1.2) under the same 
discounted rate (0% per month), genetic superiority in terms 
of HL due to Ig increased from 45.4 to 47.9 days 
(approximately 5%). This result indicates that increasing the 
discounted rate slightly shortens HL whereas increasing the 
survival rate greatly prolongs HL. 
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Applying SNP to maximize the economic weights of later 
lactation milk yields 

The net merit, selection accuracy, and genetic 
superiorities associated for selected bulls as breeding 
candidates based on selection index (Ig) are shown in Table 
4, for which SNP was applied to increase the reliabilities of 
the trait- specific EBVs. Net merit and selection accuracy 
decreased with increasing survival rate with or without 
input from SNP. In the first scenario we tested, the ratio of 
net merit at 20% increased survival rate to that of the 
current survival rate was 0.9823 to 0.9830 without applying 
SNP but 0.9863 to 0.986 with SNP. In particular, the net 
merit ratio based on the second scenario, in which the 
reliabilities of second- and third-lactation traits and HL 
were increased by 20% while those of first-lactation traits 
were maintained at the current levels, was relatively high, 
i.e., 0.9890 to 0.9897. Similarly, the ratio of selection 
accuracy at 20% increased survival rate to that of the 
current survival rate was 0.9801 to 0.9802 without SNP but 
0.9833 to 0.9845 with SNP, with that of the ratio based on 
the second scenario being particularly high, i.e., 0.9861 to 
0.9873. Therefore, using SNP reduced the decreases in net 
merit and selection accuracy that otherwise were 



accompanied by the decrease in reliabilities of later 
lactation traits. In all three scenarios we tested, appUcation 
of SNP regarding later lactation traits and HL increased net 
merit and selection accuracy. This trend becomes more 
evident as the difference in reUabiUty between second- and 
third-lactation traits and first lactation traits becomes larger 
(as in the second scenario we tested). 

In all alternative scenarios tested, applying SNP 
increased the first five lactation milk yields compared with 
these traits in the current scenario without SNP, and this 
increase became more pronounced as reliability increased 
(third scenario). Similarly, adding SNP data increased HL in 
all the scenarios (particularly the second) compared with 
that of the current scenario without SNP. The longer HL that 
occurred during the second scenario compared with that 
which occurred during the third scenario, where the 
reliabilities of second- and third-lactation traits and HL 
were increased by 20% while those of first-lactation traits 
were increased by 10%, reflects the genetic correlation 
between HL and the first-, second-, and third-lactation milk 
yields (-0.006, 0.118, and 0.257, respectively). Therefore, 
applying SNP to increase the reliabilities of EB V of second- 
and third-lactation traits and HL while maintaining those of 



Table 4. Net merit, selection accuracy, and genetic superiorities of selected bulls as breeding candidates for first through fifth lactation 
milk yields, and for herd life based on selection index (Ig)' at increased survival and discount rates when reliability of EBV is increased 
by SNP^ 





Discounted rate 


Survival rate' 


l6 


1.1 or 1.2 times reliability by SNP 




(%/month) 


(without SNP) 


1 to 1.1" 


1 to 1.2'* 


1.1 to 1.2* 


Net merit' 


0/12 


1 


663.5 


687.8 


710.0 


719.2 




0/12 


1.2 


651.8 


678.4 


702.2 


709.2 




0/12 


Ratio of 1.2/1 


0.9823 


0.9863 


0.9890 


0.9861 




4/12 


1 


665.9 


689.7 


711.4 


721.2 




4/12 


1.2 


654.6 


680.6 


704.1 


711.6 




4/12 


Ratio of 1.2/1 


0.9830 


0.9868 


0.9897 


0.9867 


Selection accuracy 


0/12 


1 


0.806 


0.836 


0.863 


0.874 




0/12 


1.2 


0.790 


0.822 


0.851 


0.860 




0/12 


Ratio of 1.2/1 


0.9801 


0.9833 


0.9861 


0.9840 




4/12 


1 


0.810 


0.839 


0.865 


0.877 




4/12 


1.2 


0.794 


0.826 


0.854 


0.863 




4/12 


Ratio of 1.2/1 


0.9802 


0.9845 


0.9873 


0.9840 


A total of the 1st to 5th 


0/12 


1 


2,717.9 


2,830.5 


2,925.3 


2,959.4 


lactation milks (kg) 


0/12 


1.2 


2,699.3 


2,812.0 


2,906.8 


2,940.0 




4/12 


1 


2,721.6 


2,834.2 


2,929.0 


2,963.3 




4/12 


1.2 


2,704.0 


2,816.6 


2,911.4 


2,944.9 


Herd life 


0/12 


1 


45.4 


48.3 


51.0 


50.3 


(AHL,d) 


0/12 


1.2 


47.9 


50.7 


53.5 


52.9 




4/12 


1 


44.9 


47.7 


50.5 


49.8 




4/12 


1.2 


47.3 


50.1 


52.9 


52.3 



' Is = biMi+b2M2+b3M3+b4Pi+b5P2+b6P3 where Mi and Pi are i* lactation milk yield and persistency EBV, respectively, and hi and hi is EBV for herd life 
and selection index coefficient for i"" component trait, respectively. 
^ SNP = Single nucleotide polymorphism. ' t times as high as the current survival rate. 

" Reliability increase as t times greater than the current reliability of estimated breeding value for the 1st, 2nd and 3rd lactation traits, respectively. 
' Net merit is a linear combination of the first five-lactation milk yields and herd life weighted by 1 and 0.389 in unit of genetic standard deviation. 



Togashi et al. (2012) Asian-Aust. J.Anim. Sci. 25:1073-1082 



1081 



first-lactation traits at their current levels increased net 
merit and HL by maximizing the increased economic 
importance of later lactation traits resulting from increased 
survival rates, by reducing the effect of first-lactation milk 
yield on HL, and by augmenting the effects of the second- 
and third-lactation milk yields on HL. 

CONCLUSION 

Net merit, selection accuracy, and HL increased by 
adding lactation persistency to the selection index whose 
omponent traits were only milk yields. Lactation 
persistency of the second and (especially) third parities 
contributed to increasing HL while maintaining the first five 
lactation milk yields compared with the selection index 
whose only component traits were milk yields. A selection 
index comprising the first three lactation milk yields and 
persistency accounted for 99.4% of the net merit derived 
from a selection index whose components were identical to 
those for net merit. We consider this selection index to be a 
practical means to increase lifetime milk yield in the 
absence of data on HL. 

Increasing survival rate increased the first five lactation 
milk yields during the 120 months from the start of the 
breeding scheme. In particular, increasing survival rate by 
20% increased the total discounted milk yields of the first 
through fifth lactations during the 120 months by 1.57 times 
compared with the total discounted milk yields at the 
current survival rate. It is important to improve survival rate 
through both genetics and management. 

Survival rate increased the relative economic 
importance of later lactation traits. However, the reliabilities 
of later lactation traits are lower than those of earlier 
lactation traits, and applying SNP to improve the 
reliabilities of the EBVs of later lactation traits and HL 
made the most of the increased relative economic 
importance of later lactation traits that resulted from 
increases in survival rate. Overall, using SNP to increase the 
reliabilities of EBVs of second- and third-lactation traits 
while keeping those of first-lactation traits at their current 
levels increased net merit and HLby making the most of the 
increased economic importance of later lactation traits 
subsequent to increase in survival rate, by reducing the 
effect of first-lactation milk yield on HL {rQ= -0.006), and 
by augmenting the effects of second- and third-lactation 
milk yields on HL {tq = 0.118 and 0.257, respectively). 
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